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SUMMARY 


This theoretical study is carried out to investigate the effect of radiation on ice crystal 
spectrum in the upper troposphere. First, an explicit expression is obtained for the ice 
crystal growth rate that takes account of radiative and kinetic effects. Second, the 
expression is used to quantitatively analyze how radiation broadens the ice crystal 
spectrum and then reveal a new precipitation mechanism in the upper troposphere and the 
stratosphere. Third, the radiative effect is used to explain the subvisual clouds near the 
tropopause. 
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1. Introduction 


Tropical precipitation is still a challenge in meteorology (Simpson et al. 1988). It 
involves multi-scale processes. Small-scale precipitating convection has a propensity to 
spontaneously organize into coherent structures or cloud clusters (e.g., Ludlam 1980; 
Moncrieff 2004), and cloud clusters in turn generate cloud shields in their mature and 
later stages (Zipser 1969; Houze 1977, 1982). Since cloud shields are widespread (or 10 4 - 
10 > km 2 in area; Houze 1982; Nesbitt et al. 2000), they modulate atmospheric radiation 
and thus change large-scale vertical circulations (Houze 1982). As a feedback, large-scale 
circulations modify clouds and precipitation (Raymond and Zeng 2005). In this multi- 
scale feedback, convective clouds and large-scale circulations are highly coupled through 
upper-tropospheric cirrus (e.g., Raymond and Zeng 2000). 

Theoretically speaking, cloud-resolving models can be used to evaluate the cirrus in 
the coupling between convective clouds and large-scale circulations. However, recent 
evaluations of cloud-resolving model simulations showed that the models predicted 
excessive clouds especially in the upper troposphere (e.g., Petch and Gray 2001; Zeng et 
al. 2007). Moreover, Stan - and Quante (2002) compared the results of cirrus models using 
bin and bulk representation of microphysical processes. They pointed out that the bulk- 
parameterization model diagnosed a significant population of large ice particles while the 
bin model was dominated by small ice particles. Those modeling studies suggest that 
some physical processes be overlooked in the current cloud models. 

Recent aircraft observations of ice particles (e.g., Heymsfield 1986; Heymsfield et al. 
2006) have provided clues as to the processes that were overlooked in the current cloud 
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models. Heymsfield (1986) showed that long-lasting subvisual clouds near the tropical 
tropopause consisted of ice crystals with small size and low number density. Heymsfield 
et al. (2006) reported that the maximum measured ice particle diameter increased with air 
temperature from -70 to -85°C when convection existed below but decreased in the 
absence of deep convection. 

Those modeling and observational issues may involve cloud-radiation interaction (see 
section 6 for discussion). Previous studies have proposed that radiation affects the 
diffusional growth of ice crystals (Hall and Pruppacher 1976; Stephens 1983; Hallett 
1987; Hallett et al. 2002). Stephens (1983) showed that radiation significantly increases 
the survival distance of falling ice crystals in an unsaturated environment. He further 
suggested that the radiative effect on ice crystal growth be introduced in cloud models. 
However, few current cloud models (Wu et al. 2000) account for it. One possible reason 
is that the ice crystal growth rate with the effect is expressed implicitly through the ice 
crystal surface temperature, making it less easy for a cloud model to adopt. Hence, it is 
imperative to obtain an explicit expression for the ice crystal growth rate first. 

Accordingly, the radiative effect is studied in the present paper beginning with an 
explicit expression for the ice crystal growth rate. The expression is then used to 
quantitatively analyze how radiation broadens the ice crystal spectrum and finally 
propose a new precipitation mechanism in the upper troposphere. The paper consists of 
seven sections. Section 2 describes a theoretical framework for ice crystal growth that 
takes account of radiative and kinetic effects. Sections 3 and 4 introduce linear and 
nonlinear models to show how radiation affects ice crystal spectrum, respectively. 
Section 5 presents a condition for haze aloft. Section 6 uses the present results to explain 
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some clouds in the upper troposphere. Section 7 concludes. Except when specified, the 
paper follows the symbol definition in appendix A. 

2. Theoretical framework 

2.1 Heat balance of an ice crystal 

The heat balance of an ice crystal in a rarefied atmosphere is represented with the 
following form (Roach 1976; Hall and Pruppacher 1976; Stephens 1983; Hallett 1987; 
Hallett et al. 2002): 

L s AnCDFp \p„ (T s ) - p v ]f m = 4 nCKF a ( T e - T s )f g + SF s - Ss 0 a(T s 4 - e r r,T, 4 ) , (2.1) 

where T s and T e refer to the ice crystal surface and the environmental air temperature, 
respectively; p and p si denote the density of the ambient water vapor and saturated water 
vapor near a crystal, respectively; S denotes crystal surface area; f„ and /q are the 
ventilation factors expressed in Hall and Pruppacher (1976). Other factors are: 

(1) C - the stationary diffusion shape factor that accounts for the influence of crystal 
shape on the water vapor field around a crystal. For a spherical crystal of radius r, C=r. 

(2) D - the diffusivity of water vapor in air that is described by (Hall and Pruppacher 
1976) 

D = 2. 1 1 X 10~ 5 (r/r 0 )■ 9 V 0 //?) [mV 1 ]. (2.2) 

(3) K - the coefficient of thermal conductivity of dry air that is calculated with the 
formula of Beard and Pruppacher (1971) 

K = [2.38 + 0.00703(7 - 7 0 )]x 10" 2 [J m’ 1 s' 1 K 1 ] , 

which is independent of air pressure even when air pressure is low (N. Fukuta 2004, 
personal communication). 
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(4) Fp — a factor that accounts for the kinetic effect on water vapor diffusion. It is 


expressed as (Hall and Pruppacher 1976) 

F„=r*/(r*+ O 


where 


r =S/4nC 


r = 


f ~ V /2 

2 n 




Qfm 


2p{i-py x 


(2.3) 

(2.4) 

(2.5) 


(5) F a - a factor that accounts for kinetic effect on heat diffusion. It takes the form 
(Fukuta and Walter 1970; Carstens 1972; Carstens et al. 1974; Hall and Pruppacher 1976; 
Pruppacher and Klett 1997) 


F.=r/{r+f Q ), 


( 2 . 6 ) 


where 


4= 


( „ V/2 

In 


kRJ,j 


Kf 0 


p a C p 2a{2-ccy 


(2.7) 


The factors 2(3(2-/3) A and 2a{2-a)' x in (2.5) and (2.7) are referred to as the apparent 
deposition and thermal accommodation coefficients, respectively (Chapman and Cowling 
1970; Fukuta 2004 1 ). Based on experimental measurements (e.g., Fukuta and Walter 
1970; Shaw and Lamb 1999), the thermal accommodation coefficient «=0.7, and the 


1 In a copy of his manuscript “Advanced Cloud Physics” received in 2004, Prof. N. 
Fukuta used the representative deposition coefficient and thermal accommodation 
coefficient to name the two factors, respectively. 
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deposition coefficient J3 = 0.04exp[-(7- T 0 )/85] that fits the data of Fukuta and Walter 
(1970) when 7>-85°C. 


(6) F s - the solar flux absorbed by an ice crystal. 

(7) £<3 — the bulk absorption efficiency of an ice crystal for blackbody radiation. It is 
defined as the ratio of the infrared flux absorbed by a crystal to the radiative flux emitted 
by a blackbody at the crystal surface temperature. Based on the results of Stephens 
(1983), the efficiency is approximated with s 0 =0.9[l-exp(— r*/6)] for a crystal surface 

temperature of T re f= 26 IK, where r is measured in micrometers. The Wien displacement 
law shows that the wavelength of peak emission for a blackbody is proportional to the 
reciprocal of temperature. Since radiation absorption depends on the dimensionless size 
parameter 2nr IX where X is radiation wavelength, the bulk absorption efficiency at 
temperature T is approximated with 

e 0 (r\T) = 0.9[l-exp(-r//6r re/ )]. (2.8) 

(8) s r -the relative bulk absoiption efficiency of an ice crystal for infrared radiation. It 
is defined as dsa where £ represents the ratio of the infrared flux absorbed by a crystal to 
the infrared flux incident on the crystal. Let T m j„ and T max denote the minimum and 
maximum temperature in the atmosphere and its underlying surface. The radiation 
emitted by a blackbody at temperature T n ,j„ and T max may impinge on a crystal as two 
extreme limits. Hence, 


£ o . ( r < g (r * J) < 

£ 0 (r,T) 


gpC^^max) 
So(r\T) ’ 


which becomes 

T nnn/ T < S r < ^ma J T 
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with the aid of (2.8). Obviously, s r varies around one. For the sake of simplicity, s r = 1 is 
assumed except in section 5. 

(9) ;/ -the ratio of the infrared flux incident on crystal surface to the radiative flux 
emitted by a blackbody at the environmental air temperature. It is approximated with 


n= 


F + + T 


(2.9) 


2 < 9 

where F~ and F represent the upward and downward fluxes of infrared radiation, 
respectively. Given a vertical distribution of atmospheric variables, the fluxes and thus 77 
are determined with radiative processes (e.g., Fu and Liou 1992; Chou etal. 1995). 


2.2 An explicit expression for ice ciystal growth rate 

The growth rate of an ice crystal is expressed in terms of ice crystal surface 
temperature T s (e.g., Hall and Pruppacher 1976). That is 

^ = 4 sCDF, [A. - PATM» . (2-10) 

where m is crystal mass. Next, (2.1) is used to eliminate T s in (2.10) to derive an explicit 
expression for dm/dt. 

To simplify the explicit expression and its derivation, a critical case of an ice crystal in 
a specific environment is considered first. In the environment, air temperature is still T e , 
but humidity is so specific that the crystal does not grow or shrink. Let the symbol T sc 
denote the crystal surface temperature, and e c and Hj C denote the partial water vapor 
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pressure and relative humidity in the critical case - , respectively. Thus, the crystal in the 
critical case has a heat balance between radiative cooling and conductive heating. That is 

& 0 °(T; C -e^:)-SF s = 4nCKF a (T e -TJf Q . (2.11) 


Since water vapor is spatially uniform in the critical case, 

PjT sc ) = e c (RX) 1 - 

After introducing the symbol 

7 = T, c )/T e 


and the approximation 


(2.1 1) is solved with 


T s >T e 4 +4T e \T sc -T), 


( 2 . 12 ) 


(2.13) 


(2.14) 


r = 


x 


where the symbol" 


Z = 


4(1 + /) 


st^r; 


1 - e r 7j - 


£ o°r e j 


(2.15) 


(2.16) 


*ckfj q 

Using (2.12) and the ideal gas law for water vapor, the critical relative humidity 
Hj C =e c /E S i(T e ) is rewritten as 


The critical relative humidity H ic is different from the effective radiative relative 
humidity defined by Hall and Pruppacher (1976). The latter is associated not only with 
the radiative effect on ice crystal growth but also the diffusivity of water vapor in air. 

The symbol /, as introduced in (2.20), represents the ratio of the infrared flux to the 
conductive heat flux from air to an ice crystal, where both air and ice crystal are treated 
as black bodies. 
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(2.17) 



using (2.14), (2.16) and the ideal gas law for water vapor. Equation (2.20) is still 
equivalent to (2.1) because the known variable T sc always satisfies (2. 11). 

Fortunately, (2.20) is formally equivalent to the energy equation of the same crystal in 
an environment that is characterized as having an air temperature T sc , partial water vapor 
pressure (eT sc /T e ), coefficient of thermal conductivity K{ 1+/), and no radiative effect. 
Deriving the explicit expression for ice crystal growth rate [or the left side of (2.20)] just 
as done in many textbooks on cloud physics (e.g., Rogers and Yau 1989) gives 

^ = [(fZkZZi!—] 1 (2.2 1) 

dt A(T c ) + B(T c )l EJTJ J 

where 
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( 2 . 22 ) 


A(T) = 


L. 




FJ 0 K(\ + X )T R r T 


B(T) = 


RJ 

Fpf m DE si (T) 


(2.23) 


To remove T sc from (2.21), A(T SC ) and B(T SC ) are changed to 


A(T sc )*A(T e )^ 


II 

t: 


B{T sc )*B{T e ) 


T e E sl (T e ) 
T sc E si (T sc ) 


with the aid of (2.2) and \T sc -T e \«T e . Substituting (2.17) and the preceding two 

approximations into (2.21) yields an explicit expression for the growth rate of an ice 

crystal in a general environment. That is 

dm _ 4nC(H : - H k ) 
dt (l-r) 2 [H lc A(T e ) + B(T e )]’ 


which does not involve T s . Expression (2.24) shows that an ice crystal grows when 
relative humidity Hi>H ic and shrinks when //,<//, c . If £ r rj= 1 and F s = 0, H ic = 100%, which 
reduces the expression into the traditional one without the radiative effect (e.g., Rogers 
and Yau 1989). 

Infrared and solar radiation effect H ic . Due to y«l and %«\, (2.15) and (2. 18) show 
that 


(1 ~H lc )cc 


\~E r T]- 


e^F £ j 


Since infrared absorption of ice crystals is more important than solar absorption, the 
effect of solar radiation on ice crystal growth is neglected in the paper (i.e., 
fJ£ 0 oT* =0). Next, Hic is discussed by dividing (1 into two factors: (l-^, ;;) and 
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(1 -H ic )l(\-£ r Tj). The former factor involves radiative processes and the latter one depends 
on crystal size, shape and other non-radiative variables. 

2.3 Infrared ratio ij 

The infrared ratio ;; approximately measures the relative magnitude of inward and 
outward infrared fluxes on the surface of an ice crystal. When £ r 7j<\, H ic < 100% and thus 
radiation enhances crystal growth. Otherwise, H ic > 100% and radiation benefits crystal 
shrinkage. Since /; varies with the variables on air temperature, water vapor and clouds, 
the following example is presented to show how ;; varies with height through the 
variables. 

Consider the tropical atmosphere observed during the Kwajalein Experiment 
(KWAJEX) that was conducted around Kwajalein Atoll from July 23 through September 
15, 1999. Figure 1 shows the mean air temperature and mixing ratio of water vapor 
versus pressure. Air temperature decreases considerably with pressure from 300 to 100 
hpa. 

Given a cloud distribution in the atmosphere, /; is calculated using (2.9) and the code 
for radiative transfer (Chou et al. 1995). Its results are displayed against pressure in Fig. 2 
for three given cloud distributions. When the sky is clear, ij >1 above the 250 hpa level 
and // <1 below that height. The factor 77 reaches a minimum of 0.89 at 428 hpa and a 
maximum of 1.8 at 100 hpa. Comparing Figs. 1 and 2 shows that the height for the 
maximum 7; corresponds to the height for the minimum temperature. Hence, 77 is largest 
at the tropopause. 


11 



Consider a thick stratiform cloud with the top at 600 hpa. Since the cloud behaves like 
a black body, /; is close to one in the cloud. Above the cloud, ?j varies with height similar 
to in a clear sky. Nevertheless, 77 is smaller than that for a clear sky, because the upward 
infrared flux decreases with decreasing the underlying-blackbody temperature from the 
sea surface temperature to cloud top temperature. 

Suppose the cloud top increases from 600 to 200 hpa. The cloud is so thick that it 
blocks the upward infrared radiation from the sea surface and the air below. Thus, it 
radiates upward as a black body but at a relatively low cloud top temperature. Hence, 
high cirrus significantly decreases 77 and therefore greatly enhances the growth of ice 
crystals above itself. 

2.4 The effect of crystal size and shape on H ic 

Once Gyij deviates from one, the radiative effect on ice crystal growth is proportional 
to a ratio of (\-Hi C )/(l-£ r Tj)- The ratio depends on crystal size and shape as well as air 
temperature and pressure, which is discussed next. 

Consider a spherical crystal. Given T and p, varies with radius. Figure 

3 displays the ratio versus radius when fg= =1. The ratio increases with crystal size. Hence, 
Hfc decreases and increases with crystal size when rj< 1 and rf> 1 , respectively, because s r 
is close to one (see section 5 for the effect of s r on H ic ). 

The ratio (\-H, c )/(l-£ r r/) also varies with height through air temperature and pressure. 
As shown in Fig. 3, the ratio increases considerably with decreasing height. It, for 
example, increases from 59 to 81% at ;• =1000 pm when (7, p) changes from (-80°C, 100 
hpa) to (- 6 °C, 500 hpa). 
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Figure 3 is suitable for non-spherical crystals, too. A columnar-shaped crystal is 
modeled by a prolate spheroid of semi-major and minor axis lengths a and b, for which 
(Pruppacher and Klett 1997) 


C = Vo 2 -ft 2 /in [(a + Ja 2 -b 2 )/b) . (2.25) 

A plate-like ice crystal is modeled by an oblate spheroid of semi-major and minor axis 
lengths a and b , for which (Pruppacher and Klett 1997) 

C = Vtf 2 -b 2 1 / arcsin y/l-b 2 la 1 . (2.26) 

After defining r* in (2.4) and using it as the horizontal axis. Fig. 3 displays as H ic versus 
r for non-spherical crystals. 


3. Two unsteady modes 

3.1 A linear model 

Expression (2.24) is used to analyze how radiation impacts ice crystal spectrum for a 
given relative humidity. Consider an ice crystal that deviates about its initial state; then 
(2.24) is linearized as 

ds/ dt = k(s 0 )[H l - H ic (s 0 )] + {s - s Q )/r (s 0 ) , (3.1) 


where subscript “0” indicates the initial state and 
s = 4/Z7 - * 2 

T(s) = -[k(s)dH ic /ds]~' 

k(s) = A-i 

V-/) 2 [H lc A(T e ) + B(T e )\ ds 


(3.2) 

(3.3) 

(3.4) 


Buoyancy waves and underlying clouds modulate Hj and H IC , respectively. Without 
loss of generality, suppose that 
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H i - H ic( s o) = + A H sin («*) 


( 3 . 5 ) 


where the forcing parameters AHto, Ah and co are constant. Substituting (3.5) into (3.1) 
and then solving the resulting equation yields 


kA, 


S — Sr, = “ 


T 2 + C0 2 


\t 1 sin(ftf) + cD(cos(rot) - 1)] 


\ + {cot) 


(3.6) 


where k and r take their values at s=So. As shown in the preceding equation, crystal 
growth accelerates when A H > 0 and A// /0 >0, and shrinkage is enhanced when Ah < 0 and 
AHm<0. These two unsteady modes broaden and narrow the ice crystal spectrum when 
rj< 1 and /;> 1, respectively, which is detailed in section 4.2. 

The third term on the right hand side of (3.6) represents the bulk ice crystal growth in 
a period due to a buoyancy wave or underlying-cloud variation. It has a factor of 
{cot)/[\ + (cot) 2 \, which reaches a maximum at ov=\. For a given r, the factor (or the 
spectrum broadening) is directly proportional to the forcing period 2n/co when cot» 1 and 
inversely proportional to the forcing period when coz« 1 . 


3.2 Timescale for spectrum broadening 

The second and third terms on the right hand side of (3.6) share a common factor of 
r(e ,/r - 1) that is inversely proportional to r. For solid spherical crystals without involving 

radiation, r— »oo, which results in r(e' /r — 1)=0. Hence, an ice crystal spectrum with |zj<oo 

is broadened or narrowed with respect to a spectrum with r— »co. In brief, r measures the 
spectrum broadening (or narrowing) due to radiation. 
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The timescale r relies on s r , ij , crystal size, shape, air temperature and pressure. As 
shown by (2.15), (2.18) and (3.3), r is almost inversely proportional to 1 -£ r ?j. Thus, the 
next discussions address r(l-£ r rj) versus crystal size and shape first and then air 
temperature and pressure. 

Figure 4 displays z{\-£ r rj) versus r* for five crystal shapes when the air temperature 
and pressure equal -80°C and 100 hpa, respectively. Using ice crystal parameters from 
observations, the timescale is calculated with (3.3), which is detailed next. 


• A graupel particle is treated as a sphere with a bulk ice density of 0.9 g kg' 1 and 
ventilation factors f„ and fg are assumed to be ones. This crystal shape is so simple that r 
versus r is easily understood with (3.2)-(3.4). Thus, the shape is discussed here as a 
reference for other shapes. Figure 4 displays the timescale T{\-£ r ij) of graupel against r*. 
The timescale decreases and increases with r when r is small and large, respectively. It 
reaches a minimum of 3.5 days at r*= 5 pm. 


• A bullet-shaped crystal is modeled with a prolate spheroid of semi-major and minor 
axis lengths a and b , for which the surface area 

. ^a 2 -b 2 


S = 2ttH 


b + 


V a 2 -b A 


arcsin- 


(3.7) 


The expressions for the mass and terminal velocity of bullets come from Heymsfield and 
Iaquinta (2000); the expression for the dimensional relationship from Heymsfield (1972); 
and the expressions for the ventilation factors from Hall and Pruppacher (1976). Using 
those expressions, the timescale for single bullets is calculated and shown against r in 
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Fig. 4. The curves of timescale versus r are not smooth because the curves for crystal 
parameters (e.g., mass, dimensional relationship) versus size are not smooth either. 
Generally speaking, the timescale curve for single bullets is similar to that for graupel 
except that z(\-£ r Tj) reaches a minimum of 3.3 days at r*= 12 pm. 


• A plate-like crystal is modeled with an oblate spheroid of semi-major and minor axis 
lengths a and b, for which 


S = 2nd 


a + 


Vtf 2 -b‘ 


In 


Vo 2 -6 ; 


+ a 


(3.8) 


The mass of the plates with dendritic extensions is calculated with the bulk density of 
Heymsfield (1972). The dimensional relationship of the plates is determined with the 
observations of Auer and Veal (1970). The tenninal velocity is calculated with the 
expression of Heymsfield (1972). With those relationships, the timescale for plate-like 
crystals is calculated and shown against r* in Fig. 4. In contrast to graupel, plate-like 
crystals have a similar timescale versus r except that z{\-£,if) reaches a minimum of 2.5 
days at r = 8.9 pm when /-*<200 pm. The decrease in timescale at /-*=200 pm coincides 
with the emergence of dendritic extensions at that size. 


• A column-like crystal is modeled with a prolate spheroid. The expressions for the mass 
and tenninal velocity of columns come from Heymsfield and Iaquinta (2000), and that for 
the dimensional relationship from Auer and Veal (1970) and Heymsfield (1972). Figure 4 
displays the timescale for column-like crystals against r and shows that r(l-£;.//) is close 
to 3.6 days when r is between 10 to 60 pm. 
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• A rosette crystal consists of bullets. It is modeled approximately with a sphere. In the 
present paper, rosettes with four bullets are assumed unless specified. The mass and 
terminal velocity of rosettes are calculated with the formulas of Heymsfield and Iaquinta 
(2000). Figure 4 also displays the timescale for the rosettes against r and shows that r( 1- 
s r rj) reaches a minimum of 0.75 day at r =44.8 pm. 

The timescale for spectrum broadening, as summarized in Fig. 4, varies with crystal 
shape. It varies with crystal size similarly for different shapes. Moreover, it varies 
considerably with height via air temperature and pressure. Figure 5 displays the timescale 
versus r when (T, p) equals (-42°C, 250 hpa) and (-6°C, 500 hpa), respectively. In 
comparison to Fig. 4, the figure shows that the timescale versus r* is similar at different 
heights, but the timescale increases significantly with height. These characteristics on the 
timescale are helpful in understanding the numerical results in the next section. 

4. Numerical simulations 

The preceding linear model deals with the radiative effect while Hj is given. Since 
crystal growth consumes water vapor, Hj changes with time. Next, a nonlinear model that 
uses Hj as a prognostic variable is introduced to show the radiative effect in a general 
framework. 

4.1 A nonlinear model 
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Consider an air parcel that moves vertically. Based on the conservation of energy and 


water, its relative humidity is governed with 


dH _ dm 


dt 


dt 


(4.1) 


where m, (j= 1 , 2, ...) denotes the mass of different sized ice crystals and 


am=- 


L,g g 


l KC p T R d 


(4.2) 


Q 2 cn= Pa 


RJ , R d L, 


+ 


EAT) pRs-CpTj 


(4.3) 


The expressions for Q\ and Qi are derived similar to that in Rogers and Yau (1989, pl06) 
for water vapor condensation. Since the present crystal growth rate takes account of the 
radiative effect, this parcel model represents the radiative effect on ice crystal spectrum. 
To model the spectrum evolution accurately, (2.24) is changed to 

da 2 _ AnC(H l - H :c ) dm . 
dt (\-y) 2 [H, c A(T) + B(T)Yda 2> ' 

which uses a 2 rather than a as a prognostic variable because the tenn on the right hand 
side of (4.4) is almost constant for spherical crystals. In the following numerical 
simulations, (4.1) and (4.4) are used as prognostic equations and are integrated with 1024 
bins. 


4.2 A mechanism for precipitation formation 

Consider a stationary air parcel with s r = 1 and a constant Based on the linear model 
in the preceding section, the evolution of the ice crystal spectrum is reasoned as follows. 
Since (\-Hj C )l(\-rj) increases with crystal size, Hi changes so that eventually it lies 
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between H ic for small and large crystals, no matter what the initial relative humidity is. 
When tj< 1, Hj is lower than H ic for small crystals but higher than that for large ones. As a 
result, small crystals shrink due to sublimation with some eventually disappearing. Large 
crystals grow due to deposition with some of them becoming precipitating crystals. This 
sublimation-deposition process broadens the ice crystal spectrum, and therefore converts 
cloud ice to precipitation. In contrast, when ij> 1, Hj is higher than H ic for small crystals 
but lower than that for large ones. As a result, large crystals shrink due to sublimation 
while small ones grow due to deposition, narrowing the ice crystal spectrum. 

To support this reasoning, a default numerical experiment is done to simulate plate- 
like crystal growth in an air parcel near the tropopause. Suppose that r=-80 o C,/?=T00 hpa 
and the initial spectrum of ice crystals 

dN{a) = c x ae C2(, da (4.5) 

where N(a ) denotes the number of ice crystals with the semi-major axis length less than 
a. After choosing Ci=1.25xl0 16 m' 5 and C2=0.5 pm' 1 , the initial concentration of ice 
crystals is 5xl0 4 m" 3 and the ice water content 9x1 O' 6 g m' 3 , close to observations for 
subvisual clouds (Heymsfield 1986). Given ^=0.5 and the initial relative humidity 
Hf= 100%, the ice crystal spectrum and //, are simulated using the time step A/=8 s. The 
modeled //, is displayed against time in Fig. 6. It decreases slowly with time. Meanwhile, 
the ice water content increases slowly due to water conservation (figure omitted). 

Figure 7 shows the evolution of the mass density dM(a)/d\na from t = 0 to 20 days, 
where M(a) represents the mass of ice crystals with the semi-major axis length less than 
a. As shown in the figure, small crystals shrink due to sublimation and large ones grow 
due to deposition. Correspondingly, the median size of ice crystals increases with time. 
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Since the small and large crystals grow oppositely, Hj lies between the H ic for small and 
large crystals. Hence, Hj is close to the H jc for the median size crystals. The increase in 
median crystal size corresponds to the decrease in Hj shown in Fig. 6 . 

If no radiation is involved (or 77 = 1), the ice crystal spectrum and Hj remain at their 
initial state. Hence, radiation is responsible for the spectrum broadening shown in Fig. 7. 

To show the radiative effect when rf> 1, the default experiment is re-done after 
increasing 77 from 0.5 to 1.5 and using the final spectrum in Fig. 7 as the initial one. Since 
(4.1) and (4.4) are symmetric with respect to 77 = 1 , it is easily inferred that the spectrum 
evolution with 77 = 1.5 reverses that shown in Fig. 7. Hence, large crystals shrink due to 
sublimation and small ones grow due to deposition when ij>\. In words, radiation 
narrows the ice crystal spectrum when ij> 1 . 

4.3 Sensitivity of precipitation formation to height 

The preceding subsection reveals that precipitation can form on a timescale of days at 
100 hpa when 77<1. Here, two additional experiments with 77 = 0.5 are carried out to show 
that precipitation forms faster at lower heights for the same 77 . The first additional 
experiment uses the same parameters as the default one except for T=- 42°C and p = 250 
hpa. In the experiment, a parameter for initial ice crystal spectrum C\=\. 25x10 m’ , 
which corresponds to an initial crystal concentration of 5xl0 6 m" 3 and an ice water 
content of 9xl0" 4 g m" 3 . The experiment is performed with a time step AM). 5 s. Figure 8 
displays the evolution of the mass density from 7=0 to 6 hours, showing that precipitation 
can form on a timescale of hours at 250 hpa, much shorter than that at 100 hpa. 
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The second additional experiment uses the same parameters as the first additional one 
except for T=- 6 °C and p=500 hpa. The experiment was performed with A/=0.02 s. Figure 
9 displays the evolution of the mass density from 7=0 to 15 minutes, revealing that 
precipitation can form on a timescale of ten minutes at />=500 hpa. In summary, 
precipitation formation due to radiation is sensitive to height via air temperature and 
pressure. This sensitivity of precipitation formation to height is a result of the sensitivity 
of the timescale to height shown in Figs. 4 and 5. 

5. Spectral asymmetry between absorbed and emitted radiation 

Consider an optically-thin cloud in a clear stationary atmosphere. The cloud is so thin 
that it has almost no effect on atmospheric radiation. Thus, much of the infrared radiation 
emitted by the underlying surface can reach the crystals in the cloud through the 
atmospheric thermal infrared window, leading to an obvious spectral asymmetry between 
crystal-absorbed and emitted radiation (or Sr*l). Next, a crystal in the cloud is discussed 
as an example to show the effect of the spectral asymmetry on H ic . 

For an ice crystal in the cloud, s r is close to its maximum, or 

e r *e 0 (r\T„„)/ £o (r\T) (5.1) 

where T„ sl is the underlying surface temperature. Equations (2.8) and (5.1) show that £j> 1 
while T<T ust . 

The infrared ratio 77 , just as shown in Fig. 2, is larger than one in the upper 
troposphere but smaller than one in the middle and lower troposphere. Consider an ice 
crystal with an environment of 77 = 1 . 8 , r=-80°C, /?=100 hpa, and T ust = 300 K. Using s r in 
(5.1), Hi C is calculated and displayed against r in Fig. 10. Meanwhile, the Hj C with s r = 1 is 
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displayed in the same figure for comparison. Obviously, effects H jc considerably when 
r* <40 pm. It changes H\ c by 0.2% at r*= 12 pm. Since e r does not reverse the increase in 
H ic with increasing r*, it affects the sublimation-deposition process only quantitatively 
while rj> 1. 

In contrast, consider a crystal with an environment of 77=0.9, 7=-12°C, p=430 hpa, and 
T US f = 300 K. Using s r in (5.1), H ic is calculated and displayed in Fig. 1 1. It increases and 
decreases with r while r* is smaller and larger than i\ LiX = 1 .9 pm, respectively. In other 
words, H ic reaches a maximum of H icMAX = 1+6x 10‘ 6 % at r=r* MAX . When \<Hj<H iC MAX, 

there are two equilibrium states [or Hi=H, c (r*)\ for an ice crystal. The equilibrium state 
with the smaller r is stable that corresponds to haze aloft. The other equilibrium state 
with the larger r is unstable. A crystal at the unstable equilibrium state, once perturbed 
initially, either shrinks to reach the stable state or grows to become a precipitating crystal. 

When many ice crystals coexist under /;< 1 and some of them have H ic less than 100%, 
the sublimation-deposition process works to form precipitation just as in the atmosphere 
with £ r =\. In summary, s r cannot change the sublimation-deposition process greatly even 
while / 7 < 1 . 

6. Discussion 

The preceding sections proposed a precipitation mechanism: when /7<1, the 
sublimation-deposition process due to radiation converts cloud ice to precipitating ice 
that in turn sediments to dehumidify air. This process, different from particle 
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sedimentation, can work as a source of large particles near cloud top. Next, cloud-top 
observations are discussed to infer the existence of the process. 

6. 1 Subvisual clouds 

Heymsfield (1986) observed the microphysical structure of subvisual clouds near the 
tropical tropopause and showed that the clouds consist of ice particles with a mean 
diameter of 2.5 pm and number density of less than 5xl0 4 m°. The subvisual clouds are 
common, but not pervasive, in the Tropics and subtropics (Winker and Trepte 1998). 
They are separated from cirrus below but tend to occur with thick cirrus below 
(Heymsfield 1986; Winker and Trepte 1998). They can also persist for weeks or months 
in clear skies in the equatorial regions (Lynch and Sassen 2002). These observations 
match the characteristics of the radiative effect, which is shown next. 

In the Tropics, convective clouds drive large-scale upward motion in the troposphere 
(Riehl and Malkus 1958). Thus, horizontal water advection cannot generate subvisual 
clouds near the tropopause. On the other side, local vertical circulations (including 
turbulent motions) transport water downgradient near the tropopause. As a result, water 
accumulates near the tropopause, and thus pervasive clouds are expected there. 

The sublimation-deposition process due to radiation can break the pervasive clouds 
into subvisual clouds. Since high cirrus leads to rj<\ near the tropopause (see Fig. 2), the 
process, accompanied with high cirrus below, can convert small ice crystals to large ones 
that in turn sediment, serving as a water sink near the tropopause. Owing to the maximum 
?; and minimum T at the tropopause, the conversion of cloud ice to precipitation is the 
least efficient there. As a result, the process can separate subvisual clouds from cirrus 
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below. With the aid of the sublimation-deposition process, high cirrus can enlarge the 
crystals near the tropopause that in turn radiate upwards strongly, which resembles the 
satellite-observed positive correlation between subvisual clouds and the thick cirrus 
below. 

In the absence of high cirrus, rj> 1 near the tropopause and large crystals shrink there, 
which can explain the small size of ice crystals in subvisual clouds. Since ice crystals 
shrink on a timescale of 10 days (see Figs. 4 and 7), subvisual clouds can last for weeks 
or months in clear skies. In summary, the observations of subvisual clouds match the 
characteristics of the radiative effect, which suggests that the sublimation-deposition 
process exist near the tropopause. 

6.2 Vertical distributions of ice crystal spectrum 

Heymsfield et al. (2006) observed large crystals with a scale up to 10' pm at -85°C. 
They reported that the maximum measured ice particle diameter increased from -70 to - 
85°C when convection existed below. That observational phenomenon, involving in the 
origin of large crystals right below the tropopause, can be explained partly by the 
radiative effect. 

Suppose that convective cells eject ice crystals into the upper troposphere (Heymsfield 
et al. 2006), which are then blown by the horizontal wind to form a stratiform cloud. The 
cloud in turn imposes a perturbation on the original vertical ;; distribution. Consider, for 
example, a stratiform cloud in the atmosphere as discussed in Fig. 1. When the cloud has 
a constant cloud ice mixing ratio of 0.1 g/kg from 132 to 197 hpa, it introduces a ;; 
perturbation as shown in Fig. 12. Since ;;< 1 near cloud top, radiation can broaden the ice 
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crystal spectrum there. In contrast, rj > 1 near cloud base, which means that radiation 
narrows the spectrum there (Stephens 1983). Thus, the radiative effect contributes to the 
increase in the maximum crystal size with height. 

The stratiform cloud discussed in Fig. 12 is optically thick. It imposes a strong 
perturbation on the original vertical ij distribution so that rj < 1 in a layer. However, an 
optically thin cloud introduces a weak // perturbation so that there is no layer with rj< 1 
and therefore no layer with spectrum broadening. If convective clouds and their anvils 
arise below a thin cloud, ij above the convective clouds is decreased. As a result, ?j< 1 
near the top of the thin cloud and rj> 1 near the base, which contributes to an increase in 
the maximum ice crystal size with height and partly explains the observational 
phenomenon that the maximum measured ice particle diameter increased from -70 to - 
85°C when convection existed below. 

6.3 Cumulus cloud seeding 

Precipitation initiation in warm cumulus clouds involves the formation of large drops 
on the scale of cloud lifetimes (e.g., 15 minutes). It requires some process to accelerate 
the growth of drops at a radius of ~20 pm where neither diffusional nor collectional 
growth is significant. Its current theories (e.g., such as giant salt particles, drop collision 
due to turbulence, and cloud entrainment) still involve some uncertainties (e.g., Rogers 
and Yau 1989; Pruppacher and Klett 1997) such that the precipitation initiation is not yet 
fully understood. Based on the present study, the radiative effect is proposed as a 
candidate to initiate precipitation in warm cumulus clouds. 
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Consider a number of crystals in a blue sky. They form due to gravity waves (or other 
processes) and some of them have length longer than ~-4 pm. Since ij< 1 below the height 
of the 250 hpa level (see Fig. 2), a part of the crystals can grow to become precipitating 
ones (see section 5) that fall and then melt below the freezing level. Once the 
precipitating particles fall into cumulus clouds, they can serve as seeds to initiate frequent 
collision between drops, which would in turn generate precipitation quickly. From 
another perspective, radiation can broaden the cloud spectrum near cloud top. Following 
the same procedure with (2.24) except for variables for liquid water, an explicit 
expression for the critical relative humidity with respect to water H wc is obtained. Since 
H wc functions in a similar manner as H ic and ij< 1 near cloud top, radiation can broaden 
the cloud spectrum there and thus contribute partly to precipitation initiation. 

In the preceding three cases, the radiative effect is expected to provide large particles 
near cloud top. Based on the consistence between the characteristics of the radiative 
effect and the observations of subvisual clouds, it is inferred that the sublimation- 
deposition process exists in the upper troposphere. Using the present explicit expression 
for ice crystal growth rate, the radiative effect will be incorporated into a cloud-resolving 
model for further evaluation. 

7. Conclusions 

This theoretical study addresses the effect of radiation on ice crystal spectrum. From 
the heat balance equation of an ice crystal, an explicit expression for ice crystal growth 
rate is obtained that takes account of the radiative effect. The expression introduces three 
variables: the critical relative humidity H IC , the infrared ratio 77 , and the relative bulk 
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absorption efficiency s r . An ice crystal grows and shrinks when the relative humidity 
Hj>H ic and Hj<H ic , respectively. 

The critical relative humidity H ic is associated with 77 through an approximately linear 
relationship between (1 -H ic ) and {\-£ r if) when the effect of solar radiation is ignored. The 
ratio (\-Hic)/(l-£ r fj) depends on crystal size and shape as well as air temperature and 
pressure. Generally speaking, the ratio increases with crystal size. Thus, given s r = 1, H ic 
increases and decreases with crystal size when ;; > 1 and 77 < 1 , respectively. 

A linear model is introduced to show how radiation broadens (or narrows) an ice 
crystal spectrum, where the timescale for spectrum broadening relies on s r , ij, crystal 
size, shape, air temperature and pressure. The timescale is computed using observed 
parameters of ice crystals. It increases considerably with height through air temperature 
and pressure. Given 77 = 0 . 5 , the timescale is measured in minutes, hours and days at the 
500, 250 and 100 hpa levels, respectively. 

A nonlinear model that uses //, as a prognostic variable is introduced to simulate the 
evolution of ice crystal spectrum due to radiation. Since H ic varies with crystal size, water 
is redistributed between ice crystals with different sizes. When 77 < 1, smaller ice crystals 
shrink due to sublimation and larger ones grow due to vapor deposition, which results in 
precipitation formation even when ice water content is small. 

The infrared ratio 77 approximately measures the relative magnitude of inward and 
outward infrared fluxes on the surface of an ice crystal. Owing to the great influence of 
clouds on 77 , the precipitation formation due to radiation is associated with clouds below. 
Based on the consistence between the characteristics of the radiative effect and the 
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observations of subvisual clouds, it is inferred that the precipitation formation due to 
radiation exists in the upper troposphere. 

However, the present parcel model needs to incorporate more physical processes so 
that its results are quantitatively comparable to field observations. Introducing ice crystal 
sedimentation, for example, is important in modeling ice crystal spectrum (e.g., Jensen et 
al. 2007). A spectrally-resolved treatment of radiation (e.g., Stephens 1983; Wu et al. 
2000) is important in the accurate determination of the infrared ratio 77, the relative bulk 
absorption efficiency s r , and the absorbed solar flux F s . With the full representation of 
the physics, the model can simulate optically-thin clouds in contrast to field observations. 
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APPENDIX A List of Symbols 


a/b : semi-major/minor axis length 
C : stationary diffusion shape factor 
C p : specific heat of dry air 
D : diffusivity of water vapor 
e : partial pressure of water vapor 

e c : partial pressure of water vapor in the critical case, see (2.12) 

E si : saturation vapor pressure over ice 

/J/q : ventilation factor for mass transfer/thermal diffusion 

F s : solar flux absorbed by an ice crystal 

F* IF : upward/downward infrared flux 

FJFp : factor for the kinetic effect in heat/water vapor diffusion 

g : acceleration due to gravity 

Hj=e/E si : relative humidity with respect to ice 

H ic : critical relative humidity, see (2.18) 

K : coefficient of thermal conductivity of dry air 
L s : latent heat of sublimation 
m : mass of an ice crystal 

M(a) : mass of ice crystals with semi-major axis length shorter than a 
N{a ) : number of ice crystals with semi-major axis length shorter than a 
r =SI4nC : equivalent radius of an ice crystal, see (2.4) 

Rd/Ry : gas constant for dry air/water vapor 
p : air pressure 
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/?o=10 5 pa : reference pressure 

*2 

s=Am' : area variable for modeling, see (3.2) 

S : surface area of an ice crystal 
t : time 

T : (air) temperature 

T e : environmental temperature around an ice crystal 

T 0 : absolute temperature of ice point 

T s : surface temperature of an ice crystal 

T sc : crystal surface temperature in the critical case, see (2.1 1) 

a/p : thermal accommodation/ deposition coefficient 

y. relative change of crystal surface temperature due to radiation in the critical case, see 
(2.13) 

So : bulk absorption efficiency of an ice crystal for blackbody radiation 
Sr : relative bulk absorption efficiency of an ice crystal for infrared radiation 
cr : the Stefan-Boltzmann constant 

X : ratio of the infrared flux to the conductive heat flux from air to an ice crystal while 
both air and the crystal are treated as black bodies 
;/ : infrared ratio, see (2.9) 

t : timescale for spectral broadening due to radiation, see (3.3) 

pj p v : density of air/water vapor 

psi : saturation density of water vapor over ice 

At : time step for numerical integration 
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Caption 


Figure 1 Vertical profiles of air temperature (thick line) and water vapor mixing ratio 
(thin line) obtained from averaging 52 days of KWAJEX data. 

Figure 2 Infrared ratio j] as a function of pressure. Thick, thin, and dashed lines represent 
cases with a clear sky, middle and high clouds, respectively. 

Figure 3 Critical relative humidity H ic versus particle radius r*, where solid, thin and 
thick dashed lines represent the results when (F, p ) equals (-80°C, 100 hpa), (-42°C, 250 
hpa) and (-6°C, 500 hpa), respectively. 

Figure 4 Timescale zfl versus particle radius r for five crystal shapes when (F, p) 
equals (-80°C, 100 hpa). 

Figure 5 Same as in Fig. 4 except that (F, p) equals (-42°C, 250 hpa) and (-6°C, 500 hpa) 
in the upper and lower panels, respectively. 

Figure 6 Relative humidity //, versus time in the default experiment where F=-80°C, 
p=100 hpa, and /pO.5. 

Figure 7 Evolution of the mass density dM(a)/d\na versus the half crystal size a in the 
default experiment where F=-80°C and p= 100 hpa. The thick line denotes the initial 
spectrum, and the time interval between lines is 1 day. 

Figure 8 Same as in Fig. 7 except that F=-42°C,/?=250 hpa and the time interval between 
lines is 5 minutes. 

Figure 9 Same as in Fig. 7 except that T=-6°C, p = 500 hpa and the time interval between 
lines is 1 minute. 
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Figure 10 Critical relative humidity H ic as a function of particle radius r*, where 77=1.8, 
r=-80°C, /?=100 hpa, and T„ s f = 300 K. The thin and thick lines represent H ic when s r 
equals one and is calculated with (5.1), respectively. 

Figure 1 1 Critical relative humidity H ic as a function of particle radius r*, where ij= 0.9, 
J’=-12°C,/?=430 hpa, and T ust = 300 K. The thin and thick lines represent H ic when s r 
equals one and is calculated with (5.1), respectively. 

Figure 12 Infrared ratio ij as a function of pressure for two cases. The thick line 
represents a case for a cloud (denoted by the shaded area). The dashed line represents a 
case for a clear sky. 
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Mixing Ratio of Water Vapor (g/kg) 



Temperature (C) 

Figure 1 Vertical profiles of air temperature (thick line) and water vapor mixing ratio 
(thin line) obtained from averaging 52 days of KWAJEX data. 
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Figure 2 Infrared ratio / 7 as a function of pressure. Thick, thin, and dashed lines represent 
cases with a clear sky, middle and high clouds, respectively. 
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Particle Radius (/^m) 

Figure 3 Critical relative humidity H ic versus particle radius r , where solid, thin and 
thick dashed lines represent the results when (F, p ) equals (-80°C, 100 hpa), (-42°C, 250 
hpa) and (-6°C, 500 hpa), respectively. 
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T=-80°C, p=100 hpa 



Figure 4 Timescale T{\-s r rf) versus particle radius /• for five crystal shapes when ( T , p ) 
equals (-80°C, 100 hpa). 
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Figure 5 Same as in Fig. 4 except that ( T , p ) equals (-42°C, 250 hpa) and (-6°C, 500 hpa) 
in the upper and lower panels, respectively. 
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Figure 6 Relative humidity //, versus time in the default experiment where r=-80°C, 
p=100 hpa, and 77=0.5. 
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Figure 7 Evolution of the mass density dM(a)/d\na versus the half crystal size a in the 
default experiment where r=-80°C and /?=100 hpa. The thick line denotes the initial 
spectrum, and the time interval between lines is 1 day. 
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Figure 10 Critical relative humidity H ic as a function of particle radius r* , where 77=1.8, 
r=-80 o C, /7=1 00 hpa, and T „ st = 300 K. The thin and thick lines represent H ic when s r 
equals one and is calculated with (5.1), respectively. 
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Figure 11 Critical relative humidity H ic as a function of particle radius r*, where 77=0.9, 
r=-12°C, p = 430 hpa, and T ; „,=300 K. The thin and thick lines represent H ic when s r 
equals one and is calculated with (5.1), respectively. 
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Figure 12 Infrared ratio ij as a function of pressure for two cases. The thick line 
represents a case for a cloud (denoted by the shaded area). The dashed line represents a 
case for a clear sky. 
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Popular Summary 

Cirrus clouds impact radiation and therefore the global energy cycle greatly. However, 
their precipitation mechanism is still unclear, which leads to excessive water in the upper 
troposphere in current cloud and climate models. This theoretical study is carried out to 
investigate the effect of radiation on ice crystal spectrum in the upper troposphere. It 
reveals a new precipitation mechanism in the upper troposphere and the stratosphere, 
which well explains the subvisual clouds near the tropopause. 



